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Abstract 

Three light sterile neutrinos (z/f, vf^ and v%.) are introduced to accommo- 
date all the available neutrino data: the atmospheric neutrino anomaly is 
explained by — z/„ oscillation with maximal mixing; the solar one is due 
to u e — v s e oscillation of small angle Mikheyev-Smirnov-Wolfenstein type; the 
Liquid Scintillation Neutrino Detector data is from v e — v„ oscillations, so 
that the neutrinos can be the hot component of the dark matter. The big 
bang nucleosynthesis constraint is satisfied by taking the tau neutrino to be 
10 MeV heavy. The v T decay is discussed in a model of gauge mediated super- 
symmetry breaking. The decay mode v T — * Gj with G being the gravitino is 
proposed. The u T has a rather long lifetime ~ 10 3 — 10 13 sec. Its implication 
to the Gamma-ray Burst is discussed. 
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I. INTRODUCTION 



Various experiments have provided growing evidence that neutrinos are massive. First 
the recent data on the atmospheric neutrinos from Super-Kamiokande experiment has shown 
neutrino oscillations Q. It implies that the /x-type neutrino has a maximal mixing with 
other neutrino x (x 7^ e), and the mass squared difference is Ami. ~ 3 x eV 2 . The 
solar neutrino deficit problem can be explained by either the Mikheyev-Smirnov-Wolfenstein 
(MSW) solution |2| or the vacuum oscillation ||. The MSW solution allows two sets of 
parameters: Am^ ~ 5 x 10~ 6 eV 2 with sin 2 26 ey ~ 6 x 10~ 3 , and Am 2 y ~ 2 x 1CT 5 
eV 2 with sin 2 26„. ~ 0.8. The vacuum oscillation solution is Ami ~ 8 x 10 -11 eV 2 with 

c .y ey 

sin 26 ey ~ 0.8 Q. If the direct observation of — u e oscillations at the Liquid Scintillation 
Neutrino Detector (LSND) experiment |J is considered, relevant parameters are Am^ ~ 1 
eV 2 and sin 2 (2# e/ J ~ 10" 2 . 

In addition, astrophysics and cosmology also give us some information for neutrino 
masses. Various measurements support the existence of the dark matter. One scenario 
for the dark matter is that not all of the cosmological dark matter is cold, there is some 
hot component ||. Canonical candidate for the hot dark matter (HDM) is the neutrino 
with mass at the scale of electron Volt ||. If massive neutrinos are stable, their masses are 
limited to be less than few tens of eV from the standard cosmology, in order to avoid the 
over-closure of the Universe [[?]. 

Non-vanishing neutrino mass might be the first discovery of the physics beyond the stan- 
dard model (SM). The mass pattern of neutrinos is valuable information for the exploration 
of the physics related to the flavor puzzle in the SM. Theoretically, several mechanisms have 
been suggested to accommodate massive neutrinos. One of the most popular scenarios is the 
seesaw mechanism which naturally explains the smallness of neutrino mass by introducing 
heavy right-handed neutrinos ||. Another example is the supersymmetric extension of the 
SM with .R-parity violation. The trilinear lepton number violating interactions induce small 
neutrino masses at loop level 0. However, the large mixing of neutrinos is not expected 
from our experiences of the other SM fermions. 

Phenomenological analyses suggest the existence of light sterile neutrinos. As least one 
light sterile neutrino is necessary to explain all the three different scales of Am 2 [10]. It is 
interesting and more natural that there are three light sterile neutrinos u- (i = e, fi, r). Then 
the first two generations can fully explain all the above neutrino data. In this framework, the 
neutrino mixing pattern is uniquely fixed. The atmospheric and solar neutrino anomalies 
are due to the — u s and v e — v s e oscillations, respectively. The large mixing between 
and V s can be understood naturally. And the LSND evidence is attributed to the u e — 
mixing. At least one of the neutrino pairs should have mass around few eV, implying that 
neutrinos can be the hot component of the dark matter. This scenario was studied before 
[ TT| , [T2"| ]. In this paper we reconsider it in a different theoretical background. 



In this scheme, more physics can be addressed. The third generation neutrinos have a 
separate story. The above requirement for the first two generations leaves a wider room 



1 Recently the low mass solution has been revived which is Am^ ~ 10 7 eV 2 with large mixing 
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for the v T . The v T can be either ordinary with a mass smaller than 10 eV, or exotic, for 
instance, with m Ur ~ (1 — 10) MeV [11]]. The physics of tau-neutrinos is influenced by 
the astrophysical constraint from the big bang nucleosynthesis (BBN): the number of light 
neutrino species in thermal equilibrium at the BBN era is limited to N v < 4.2 [13|. It should 
be seriously examined in a scenario with more than four light neutrinos. It is to be noted 
that the small mixing angle for the v e — v s e oscillations causes only the active neutrinos to 
contribute at the BBN era [|TJJ]. The oscillation time is too long to make the v s e in thermal 
equilibrium. In fact, the analyses for the solar neutrino data have shown that the only viable 
pattern is the small angle MSW type Jr5| . Therefore the first two generations contribute 
a factor three to N v . As for the v T and v^, at least one of them should decouple at the 
time of BBN. If v T is light (m Ur < 10 eV), the sterile tau neutrinos must decouple from 
other particles in the primeval plasma. This is equivalent to the case of introducing just 

On the other hand, we may consider m L 



two sterile neutrinos 



> 1 MeV so that u T 



decouples [[H 
motivation 



We are interested in this latter situation which has alternative theoretical 
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In addition, the economy of the first two generation oscillation scenario 
becomes more meaningful. 

It is a theoretically challenging problem to construct a spectrum of very light sterile 
neutrinos. Several ideas have been suggested. They are related to extra symmetry [IS|, 
mirror world [IS], axino [2C], modulino [^l|], gravitino [^2], string theory ||23|| , gauge medi- 
ated supersymmetry breaking (GMSB) ||24|| , U(2) symmetry [p5| , seesaw mechanism [26| , 
compositeness f2"7fl , top-flavor model |28|] and extra dimensions p9| . 

This paper is organized as follows. In the next section, we present a mass matrix which 
is consistent with all the available neutrino data. The naturalness of the assumption is 
discussed. In Sect. Ill, the decay of the tau neutrino is studied in a supersymmetric model. 
The decay mode v T — > G7 is proposed. Further discussions including the gamma-ray burst 
(GRB) are made in Sect. IV. The summary is given in the last section. 



II. NEUTRINO MASS MATRIX 

We assume that there are three light sterile neutrinos. For each generation of fermions, 
there is an associated sterile neutrino. As will be seen in the following, the atmospheric 
and LSND data will be accommodated naturally. And the tau neutrino can be as heavy 
as (1 — 10) MeV without contradicting experiments. The full 6x6 neutrino mass matrix 
A4 consists of three parts. They are 3x3 active Majorana neutrino mass matrix M actlvc , 
3x3 Dirac mass matrix M Dirac which stands for the mixing between the active and sterile 
neutrinos, and 3x3 Majorana mass matrix of the sterile neutrinos, respectively, 

(j^active ^|"Dirac 
fl^"Dirac n^"sterile 

For the j\/f actlve , we assume that except for the v T v T component it originates from the 
ordinary seesaw mechanism ||; and the v T v T component is about (1 — 10) MeV. The ordinary 
seesaw mechanism suppresses neutrino masses to be v 2 /Mqut, where v is the electroweak 
scale and Mgut the grand unification (GUT) scale so that t> 2 /McuT ~ 10~ 3 eV. This is 
natural because there are certain evidences that the new physics beyond the SM is the GUT 
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SOfl . The introduction of the sterile neutrinos does not affect the gauge coupling unification. 
Furthermore, as can be seen below, this 10~ 3 eV mass just fits the atmospheric neutrino 
data. In other words, the neutrino mass itself is also indirect evidence of the GUT. Note 
that even without the heavy right-handed neutrinos, the seesaw mechanism still works in 
the GUT scenario. It is from the general argument based on the effective field theoretical 
point of view. The MeV tau-neutrino has been discussed in Ref. [17|j. We introduce it 
here to avoid the BBN constraint and to make the two generation oscillation scenario more 
economical. 

For the M Dirac , we adopt that it is proportional to the up-type quark masses. This simi- 
larity is understandable from the GUT point of view, if the sterile neutrinos are regarded as 
right-handed neutrinos. It is remarkable that this simple proportionality can be compatible 
with all the current neutrino data. In this quark mass matrix, the up quark is taken to be 
10~ 4 m c . The up quark mass has not yet been determined experimentally. Even a massless 
up quark is allowed JH],^]. 

For the j\/f stenle , we simply assume that it is approximately a vanishing matrix. This 
assumption makes the point of light sterile neutrinos. It could be relaxed to certain extent 
that the M stenle plays no other roles in understanding the neutrino experiments. 

Since the main understanding of the neutrinos does not involve intergenerational flavor 
mixing, without losing generality, we further assume that all these mass matrices are diag- 
onal. We make a remark here that the other zero entries of the Majorana mass matrices 
should be regarded as ~ 10~ 5 eV in general. This is inevitable from the consideration of 



quantum gravity [ 33| . However, such a number is not essential for the understanding of the 
neutrinos in this framework. 

Therefore the neutrino mass matrix Ai takes the following form, 
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(2) 



V 

where A is a constant. 

Here starts numerical analyses. As has been mentioned, v 2 /Mqut is taken to be 3 x 10~ 3 
eV. We take Xm c ~ 1 eV. It is easy to see that and v s have an almost maximal mixing with 
Am 2 ~ 3 x 10~ 3 eV 2 . This explains the atmospheric neutrino data. The mixing between 
u e and v s e is very small: sin 2 20~[2x 10~ 4 /(3 x 1CT 3 )] 2 ~4x 1CT 3 with Am 2 ~ 1CT 5 eV 2 . 
This is the small angle MSW solution for the solar neutrino data. The Dirac mass of tau 
neutrino is 100 eV. This leads to no observable physical consequences of sterile tau neutrinos 
because the Majorana mass of v T is about 10 MeV which makes the v T — v s T mixing to be 
at ~ O(10~ 5 ). As discussed before, this small mixing precludes the sterile tau neutrinos 
from contributing at the BBN era. The LSND data of the v e — oscillations is accounted 
by Am 2 ^ ~ 1 eV 2 . The intergenerational mixings of the neutrinos are determined by the 
neutrino mass matrix and the charged lepton mass matrix, which is too complicated to be 
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fixed in this framework. Since the Dirac mass of the electron neutrino is negligible, the 
v e — Up mixing depends only on M actlve and the charged lepton mass matrix. In the case 
where M actlve is diagonal in certain weak basis, we expect that the charged lepton mass 
matrix will fix the v e — mixing, namely m e /m M ~ 5 x 10~ 3 which is consistent with the 
data. The and V s have masses around 1 eV so that they can be the HDM. The v T cannot 
be stable due to cosmological considerations. 



III. THE DECAY OF v T 

In the previous section, phenomenology of three light sterile neutrinos is discussed, and 
a heavy Majorana tau neutrino is necessarily assumed. To study the 10 MeV tau neutrino 
and its decay, a theoretical model needs to be specified. First, this v T may decay via the 
ordinary weak interaction, namely into e + e~v e through the W^-boson exchange. Its lifetime 
is 

1927T 3 1 

0-3 X y- /in1 ^ m , T , r, *■ (3) 



Vt G 2 F ml\V eT \ 2 ' (m^/lOMeV) 5 )^! 2 

Since the neutrino mass matrix leads to a very small mixing, only the charged lepton mass 
matrix will determine V eT . It is generally expected that \V eT \ 2 ~ m e /m T ~ 3 x 10~ 4 which 
results in t Vt ~ 10 3 sec. However this decay mode is disfavored by the observation of the 
Supernova 1987A pi |. Even if this decay is marginally viable by assuming smaller V e 



CT ) 



something more about the v T decay can be addressed. 

There is an alternative decay channel which is irrelevant to the e — r mixing. It is 



inherent in the model of the (1 — 10) MeV v T . In Ref. [35|, a supersymmetric model for the 
charged lepton masses is suggested. In this model, because of family symmetries, r lepton 
gets mass from the nonvanishing vacuum expectation value (VEV) of the Higgs field. The 
muon gets mass only due to the sneutrino VEV which violates the family symmetries. The 
lepton number violation is allowed in the model. The sneutrino VEV, which generates a 10 
MeV mass to the tau-neutrino, is as large as (5 — 10) GeV, resulting in a mixing between 
the ordinary neutralinos and the v T WJ\. This model adopts the framework of the GMSB 



36fl . The v T in fact is the lightest neutralino (except for the gravitino) and the r the lightest 



chargino. Therefore v T can decay to G7 where the gravitino has mass of few eV's. From 



the neutralino mass matrix in Ref. [17], it can be seen that v T mixes with Bino at the level 
of m T H/M^ ~ 10~ 3 where itl t h is the mixing mass of the neutrino and neutralinos. In 
ordinary GMSB models |36| the next-to-lightest neutralino, photino 7, decays to G and 7. 



The relevant matrix element (G7I7) can be written in the same way as in current algebra 

<G 7 |7> = }(7|^?|7>, (4) 
where the supersymmetry current is 

~f = F>fG + a' u 'ZfF ia , + ... (5) 



with yF being the supersymmetry breaking scale ~ 100 TeV. From the second term in 
the above equation, the matrix element can be easily evaluated. In our model, through 
considering the suppression factor (m T ff/M^) 2 , the v T — > G7 decay rate is then given by 
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r(^G T )4— ) • c ° s ° w > ■ (6) 



This is equivalent to a lifetime of 10 13 sec, if we switch off the V eT . Even though this lifetime 
is very long compared to the BBN epoch, it is still short enough compared to the age of the 
Universe. 

The BBN constraint for the number of light neutrino species is counted by u e , and v 8 . 
A lifetime of 10 3 sec for v T is at the finishing stage of the BBN. Such tau-neutrinos may be 
considered as unstable in the BBN time. With v e being the decay product, it is probably 
helpful for the understanding of the Deuterium production and the structure formation in 
the cosmology |37|. It is however more appropriate to take v T with lifetime of > 10 3 sec as 



stable in the BBN epoch |5B |. In this case, v T decouples in the BBN time. And it must be 
Majorana neutrino f37],|38| as we have assumed. 



IV. DISCUSSIONS 



Several discussions are necessary about our scenario for neutrinos. The essential new 
ingredients of this scenario are the 10 MeV Majorana tau-neutrinos and the three light 
sterile neutrinos. A theoretical model of the 10 MeV v T and its decay have been discussed. 
However, the lightness of the sterile neutrinos is still lack of theoretical discussion. To our 
knowledge, for more than one light sterile neutrinos, there exists two attractive frameworks. 
One is the compositeness idea [27j. And the other comes from the mirror world | l9fl . They 
are worthy to be studied further. We introduce them phenomenologically in this paper. 

While the .R-parity violating supersymmetry is used to generate the MeV mass for u T , 
it would be expected that the trilinear i?-parity violating interactions contribute to the 
other entries of M actlvc . However, this contribution is rather small. It arises at one-loop 
level, where a soft term of the relevant trilinear i?-parity violating interaction, A-term, is 
involved. In the GMSB, A-terms are very small. The contribution is in general smaller than 
10~ 3 eV. 

Experimentally, our scenario can be tested in the very near future. First, the current 
laboratory mass limit for v r is 18.2 MeV p9]. It can be improved to be a few MeV in the 



near future, which could be a serious challenge to this scenario. 

Second, future neutrino experiments will determine whether the atmospheric neutrino 
deficit is due to — v s or — v T oscillation. Although the present data favor the — v T 
oscillation, the — v s „ oscillation cannot be ruled out on the basis of the global fit to the 



full set of observable PP| . Recent Super-Kamiokande experiment claims that the 



oscillation for the atmospheric neutrino anomaly is ruled out at 99% CL. [El]. However, it 



has been argued that more data and analyses are necessary before the conclusion is made 
The long baseline neutrino experiments [|3] will check the oscillation channel. On the 



other hand, the oscillation channel for the solar neutrino deficit can be hopefully picked out 



by the new SNO observation [Q together with previous experiments. 

Third, the short baseline neutrino experiment BooNE will confirm whether LSND 
anomaly is real. Even though the KARMEN [15] group has reported that a large part 
of the favored parameter region of LSND is excluded, full confirmation of the LSND results 
still awaits future experiments. For example, some of the LSND solution space is preserved 
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by a combined statistical analysis of both the LSND evidence and the KARMEN exclusion 

Cosmological and astrophysical aspects of this scenario deserve further studies. The 
implications on the BBN and HDM have been discussed in the previous sections. For the 
decay mode v T — > G7, it might be the partial reason of the GRB P?fl . In the process of the 
stellar collapse, a huge number of tau neutrinos are emitted. The decay product 7 is of few 
MeV or 100 keV which is the typical energy scale of 7 in the GRB. Suppose that the lifetime 
of v T is of order 10 13 sec ( the V er is negligible small < 10~ 7 ), 10 MeV v T emitted from the 
Supernova 1987A did not have enough time to decay when they arrived in the Earth EB 



When the supernova occurs at cosmological distance (> 10 26 — 10 27 cm) from us, however, 
the GRB will be detected. The (1 — 10) MeV v T may help our understanding of the cold 
dark matter (CDM) fl9fl . As for the sterile neutrinos, it will be harmless for our numerical 
results to assume that v s T has Majorana mass of thousand electron Volts. Then this v s T is 
also the candidate of the warm dark matter (WDM) ||50|| . 

On the other hand, if v T purely decays to G7 with a lifetime ~ 10 13 sec, too much 
cosmic background radiation and matter density are caused |7|]. For reducing the background 
radiation, one way maybe to prolong the v T lifetime from 10 13 sec until to 10 21 sec. This can 
be achieved through adjusting the supersymmetry breaking scale (along with some other 
model parameters) up to 10 4 TeV. (In this case, the gravitino mass rises from few eV to 
100 keV.) To avoid the matter density problem, we may need some non-standard cosmology. 
For example, in a recent study [^TJ, it has been shown that if the largest temperature of the 
radiation era is as low as 0.7 MeV, even a stable 10 MeV tau neutrino is consistent with the 
cosmology. 

Within the framework of GMSB, the superpartners of the sterile neutrinos can be very 
light because they are not involved in the gauge interactions. However, their mixings with 
the ordinary sneutrinos which are typically 100 GeV — 1 TeV heavy, are very small. So their 
lightness has no influence on the physics of the BBN. Meanwhile, these light scalar particles 
can be also the candidate of the CDM |52| . 

Compared to the previous studies of MeV tau neutrinos, the distinguished feature of 
this scenario is that the (1 — 10) MeV tau neutrino is Majorana particle and its decay 
products do not contain Majoron. Different from the pseudo-Dirac neutrino scenarios ||12|| , 



our explanation of the solar neutrino deficit is due to small angle MSW solution. The 
previous motivation was mainly to reduce the number of neutrino species in BBN era to be 
smaller than 3, the present one is to keep that number just as 3. 



V. SUMMARY 

Before summarizing the main points, let us list the dark matter candidates in this model. 
Besides the and V s which have been taken as the HDM, there are many candidates for the 
CDM and the WDM. They are (1 — 10) MeV z/ T 's, keV massive i^'s, light sterile sneutrinos v? 
and gravitinos. However, these are just possibilities. Because the lifetime of v T has not been 
really fixed which still ranges from (10 3 — 10 13 ) sec (or longer); the heaviness of v s T and v\ are 
just assumed, which need to be studied from the theoretical mechanism for the lightness of 
sterile neutrinos. Note that by taking the supersymmetry breaking scale as ~ 100 TeV, the 
gravitino is not the dark matter in simple GMSB models. However by raising the \fF, keV 
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G can be the WDM (regardless the R-parity violation). In addition, in the GMSB scenario, 
lightest stable baryons with mass of ~ 100 TeV in the hidden sector is also a candidate 
for the CDM ||36|| . Among the five possibilities, the long-lived and decaying v T which is 
specific in this model, is reasonably favored to be the CDM in this model according to our 
discussion. For the stable z/ T , the matter density was discussed in Ref. It is expected 
that the same discussion can be made for the very long-lived v T . The requirement Vt VT h 2 < 1 
implies m„ r > 3 MeV for the reheating temperature being 0.7 MeV. 

In summary, we have introduced three light sterile neutrinos V s and v s T to accom- 
modate all the available neutrino data. The atmospheric neutrino anomaly is explained by 

— v s oscillation with maximal mixing and the solar one by v e — v s e oscillation of small 
angle MSW type. The LSND data of the v e — oscillations with Am^ ~ 1 eV 2 can be ac- 
commodated, as providing canonical candidate of the HDM. The BBN constraint is satisfied 
by taking the tau neutrino to be 10 MeV heavy. The v T decay has been discussed in a model 
of GMSB. The decay mode v T — > G'-f has been proposed. The implication to the GRB has 
been pointed out. Although the current data do not favor this oscillation scenario, it will be 
conclusively checked by atmospheric neutrino experiments, long baseline experiments and 
solar experiments in the near future. 
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